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This work addresses nucleation rates in systems with 
strong initial recombination. Initial (or 'geminate') recombi- 
nation is a process where a dissociated structure (anion, vor- 
tex, kink etc.) recombines with its twin brother (cation, anti- 
vortex, anti-kink) generated in the same nucleation event. Ini- 
tial recombination is important if there is an asymptotically 
vanishing interaction force instead of a generic saddle-type ac- 
tivation barrier. At low temperatures, initial recombination 
strongly dominates homogeneous recombination. In a first 
part, we discuss the effect in one-, two-, and three-dimensional 
diffusion controlled systems with spherical symmetry. Since 
there is no well-defined saddle, we introduce a threshold which 
is to some extent arbitrary but which is restricted by phys- 
ically reasonable conditions. We show that the dependence 
of the nucleation rate on the specific choice of this thresh- 
old is strongest for one-dimensional systems and decreases in 
higher dimensions. We discuss also the influence of a weak 
driving force and show that the transport current is directly 
determined by the imbalance of the activation rate in the di- 
rection of the field and the rate against this direction. In a 
second part, we apply the results to the overdamped sine- 
Gordon system at equilibrium. It turns out that diffusive ini- 
tial recombination is the essential mechanism which governs 
the equilibrium kink nucleation rate. We emphasize analogies 
between the single particle problem with initial recombination 
and the multi-dimensional kink-antikink nucleation problem. 

PACS numbers: 82.60. N, 05.20, 82.30. L, 72. 20. J, ll.lOKk 



I. INTRODUCTION 

The purpose of this work is to support our approach 
to the nucleation problem of kinks and antik-kinks in the 
overdamped sine-Gordon chain at equilibrium [Q by in- 
vestigating closely related problems which are dominated 
by initial (or 'geminate') recombination processes. Ini- 
tial recombination is a process where a dissociated struc- 
ture (anion, vortex, kink etc.) recombines with its twin 
brother (cation, anti-vortex, anti-kink) generated in the 
same nucleation event [|]-f§ ■ Our theory of kink-antikink 
nucleation |IJ leads in the thermal equilibrium state of 
the sine-Gordon chain to an activation rate proportional 
to exp(— 2Eo/kT) with an activation energy of twice the 
equilibrium kink energy Eq. This can be compared with 
earlier work M which predicts a rate proportional to 
cxp(— 3Eo/kT). This significant difference in the nucle- 



ation rates is a consequence of initial recombination of 
kink antikink pairs. Although our discussion presents a 
clear physical picture, our results are apparently not ob- 
vious and have already given rise to discussions in the 
literature It is therefore necessary to further sup- 

port and explain more deeply the point of view and the 
approach taken in Ref. . 

Usually, nucleation theory jlo| (for reviews, see []Tl| , [l2|| ) 
is associated with the decay of a metastable state across 
an activation barrier. A simple picture is a Brownian 
particle which has to overcome a barrier in order to leave 
the region of attraction of a potential well (see Fig. pi). 
One of the main tasks of nucleation theory is the evalua- 
tion of the current flow out of the well for a given density 
of particles in the well. The relevant time scale of such 
an activated process is dominated by an inverse Arrhe- 
nius factor exp(Eo/kT), where Eq, k, and T denote the 
activation energy, the Boltzmann constant, and the tem- 
perature, respectively. The nucleation rate is then pro- 
portional to the Arrhenius factor. The proportionality 
constant is in the center of interest in many works and 
depends strongly on the shape of the potential and on 
the strength of the damping of the Brownian particle. In 
the present work we emphasize the particularly interest- 
ing problem of nucleation and recombination in systems 
with an asymptotically 'flat saddle' (see Fig. |l|b). We 
will consider activation out of a well onto a region where 
the potential converges fast to a constant value, without 
crossing a barrier with a maximum. Except for a few 
remarks on the underdamped limit we will restrict our- 
selves to strongly damped systems. This corresponds to 
a diffusion in configuration space on the flat saddle re- 
gion. Furthermore we will allow a very slow drift due to 
a weak external driving force. 

The main difference to the usual nucleation across a sad- 
dle with a well defined maximum consists in an enhanced 
probability of backscattering into the well. Indeed, while 
there the particle is driven away by the deterministic 
force once the barrier maximum is traversed, the par- 
ticle on a flat saddle executes a diffusive motion leading 
to a large probability of falling back into the original 
well. This process corresponds to initial recombination 
in generation-recombination kinetics and is in contrast to 
homogeneous recombination. The homogeneous recom- 
bination rate is proportional to the density of wells, and 
is much smaller than the rate of initial recombination for 
a diluted well density. The processes are illustrated in 
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the lower part of Figj^. 

In the following we mention some examples of initial 
recombination. The most prominent one appears in dis- 
sociation theory of diffusion controlled chemical reactions 
[3|. In his 1938 paper titled 'initial recombination of ions' 
[ 2 1 , Onsager determined the probability of recombina- 
tion of a pair of ions after a given initial separation r. 
It is well-known, that despite the long-range nature of 
the Coulomb force, there is a finite escape probability 
in three dimensions even at zero external driving force. 
Onsager gives for the probability of escape exp(— r a /r), 
where r Q = qiq 2 / (47re r eofcT) is called the Onsager ra- 
dius, qi,2 denote the charges of anion and cation, and 
e r £o is the dielectric permeability of the medium. While 
for small initial distances r the attractive Coulomb in- 
teraction leads to strong initial recombination, diffusion 
dominates Coulomb interaction for r larger than the On- 
sager radius which leads to a large escape probability. It 
is clear that in a dilute electrolyte and if the dissociation 
process is due to thermal nucleation, initial recombina- 
tion dominates homogeneous recombination. This has 
strong consequences on the free ion density in a super- 
imposed electric field, since initial recombination is very 
sensitive to an external force. Note that in a flat poten- 
tial, i.e., if r = 0, Onsager's result predicts a vanishing 
initial recombination rate. We will show below, that this 
is characteristic for a three-dimensional system, and it 
will be different for lower dimensions. 

A two-dimensional example is the two-dimensional 
Coulomb gas. Here, the force vanishes proportional to 
1/r which leads to a logarithmically increasing potential. 
Though the potential is not 'flat', this case is particu- 
larly interesting (see, e.g., Dissociation in the 2-d 
Coulomb gas describes, for instance, vortex-antivortex 
nucleation in superfluid Helium films. An important fea- 
ture of this system is that there is an unbinding transition 
of pairs at a critical temperature. At low temperatures 
the escape probability vanishes, whereas above a crit- 
ical temperature it is finite. This transition, which is 
known as the Kosterlitz-Thouless transition, can be seen 
as a consequence of strong initial recombination which 
inhibits free excitations below the critical temperature 
j|. We mention, that there is no simple single-particle 
picture of the Kosterlitz-Thouless transition since screen- 
ing of the excitations play a major role. 

Initial recombination in one-dimensional systems is 
faced, e.g., in photoelectric carrier generation in 1-d semi- 
conducting polymers [^). Haberkorn and Michel-Beyerle 
H considered electrons which are photo-generated close 
to an electrode in a one-dimensional conductor, includ- 
ing the image-force at the electrode. In particular, they 
discussed the relation between the total current due to a 
weak external force, and the 'nucleated' current. We will 
take up these results below, 

Note that the listed examples are all controlled by an 
interaction force rather than by diffusion alone. In this 
paper, however, we emphasize diffusion controlled, i.e., 
entropy driven, initial recombination. This requires that 



the force be short ranged. Short range forces occur, e.g., 
between kinks or domain walls in (quasi) one-dimensional 
systems. A kink-antikink interaction potential is often a 
monotonously increasing function of the separation and 
becomes exponentially flat for large separations JT3 ] . 

The paper is organized as follows. In order to fix the 
notation and for later comparison, we briefly review in 
Sect. H some standard results of the theory of nucleation 
across a barrier. In Sect. Ill we discuss the nucleation 



across flat saddles, including the effect of a weak force, 
and the initial recombination time. In Sect. IV, diffu- 



sion controlled initial recombination is illustrated for the 
example of equilibrium kink nucleation. 



II. ESCAPE ACROSS A BARRIER 

In the following we consider the motion of a particle in 
d dimensions with coordinate r, mass m, damping con- 
stant 7, and in a potential U(f). Of course, the particle 
coordinate vicariously stands for relevant variables of a 
rather large class of systems (chemical reaction variables, 
ion separation distance, liquid droplet radius, magnetiza- 
tion, etc.). Throughout this paper, we assume that the 
potential U has a minimum at r — and depends only 
on the distance r = |fj. In a one-dimensional system 
(d = 1), this means that U(r) is an even function of the 
coordinate r. The particle is furthermore coupled to a 
heat bath of temperature T. The equation of motion for 
r is a Langevin equation with a white noise force p5| . 
Equivalently, the system can be described by Kramers 
equation for the phase-space probability density /(r, v, t), 
where v is the velocity of the particle (for our purposes 
we may use the velocity instead of the momentum). Be- 
low, we consider almost always the strongly damped case 
in which Kramers equation reduces to the Smoluchowski 
equation, 



d t P + Vj = 



(1) 



which is a continuity equation for the probability den- 
sity P = J d d v / in configuration space. The probability 
current j = J d d v vf is given by 



j = (PVU + fcTVP) 

rwy 



(2) 



where the Einstein relation D — kT /rwy between diffu- 
sion constant D and damping constant 7 is related to 
local equilibrium. 

The nucleation rate J nuc can be defined as follows. 
First, one has to find a stationary solution of Kramers 
equation which is normalized in the well. Outside the 
well, this solution has to satisfy appropriate boundary 
conditions. Usually, absorbing boundary conditions are 
applied. Then, the current crossing the boundary is the 
nucleation rate. This approach has the advantage that 
the determination of the nucleation rate is reduced to the 
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solution of a stationary (time-independent) problem. 

In the remaining part of this section, we recall some 
results for one-dimensional systems and for saddles which 
are not flat in the above mentioned sense. Further- 
more, we restrict the discussion to low temperatures, i.e., 
kT <C Eq. For moderate-to-strong damping, Kramers 
found for a quadratic saddle at r — s an analytical re- 
sult which reduces in the strong damping limit to the 
Smoluchowski rate 



Jnuc = 2 LU ° L ° S exp (-E /kT) 

ZTT'f 



(3) 



The rate depends on the Arrhenius factor exp (—E /kT), 
on the damping constant 7, and on the curvatures at 
the minimum and at the saddle, luq — \J d^U(0)/m and 
ujs = \/\d^U(s)\/m, respectively. The factor 2 in Eq. 
(||) occurs because the symmetry of the one-dimensional 
potential implies two equivalent paths of escape from 
the well. In the weak damping limit (7 — > 0), Kramers 
moderate-to-strong damping result yields the transition 
state rate 



Jtst = 2 ^ exp (-E /kT) 



(4) 



This result is wrong since thermal equilibrium is assumed 
in the well which is no longer the case at very low damp- 
ing. An appropriate (slow) variable is then the action 
variable / (or the energy E). By an averaging over the 
fast angle variable, Kramers derived therefore a diffusion 
equation in action space. A solution of the nucleation 
problem in this space yields the true rate in the 7^0 
limit: 



(5) 



Here Is is the action of the separatrix trajectory at the 
saddle in phase space. We mention, that the result (||) 
remains valid also for the 'flat saddle' in Fig. 0b, since 
at 7 — > the particle escapes ballistically, i.e., without 
backscattering in the region with constant potential. 

The rates (||) and (^) are the leading order results in 
I/7 and 7, respectively. The crossover between these 
limits is known as the Kramers turnover problem. Re- 
newed attention to this problem followed Refs. [ fil||T7| 
which investigated the leading corrections to Eq. (5) and 
provided computational results illustrating the turnover. 
A particularly elegant approach to this problem was put 
forth by Melnikov |l^JT^], and independently by Risken 
et al. JT^,0 based on boundary layer theory. Later it 
has been noticed that the problem can be solved over 
the entire turnover region by normal mode analysis [^l] . 
We expect that for the 'flat saddle', the first correction 
term to Eq. (||) is strongly changed as compared to the 
quadratic saddle. In particular, there should be a strong 
dependence on the length of the diffusion region. In the 
present work, however, we consider the strongly damped 
case, or corrections to this case for finite but small I/7. 



To find the corrections, we will use a result of boundary 
layer theory. 



III. ESCAPE ACROSS A FLAT POTENTIAL 
REGION 

We consider again a symmetrical potential well. For 
decreasing r < £0, the potential is assumed to drop 
strongly and monotonously to its minimum value — Eq at 
r = 0. For r > £0, ft is assumed to increase monotonously 
and to converge fast to zero. Here £0 is a characteris- 
tic half-width of the well defined by the specific physical 
problem under consideration (see, e.g., Sect. IV). A typi- 
cal example is sketched in Fig. |l|b. As mentioned above, 
for the one-dimensional case d = I, we assume reflec- 
tion symmetry, and in higher dimensions we assume a 
rotationally invariant potential, which allows us to con- 
fine all calculations to r > 0. On a large length scale 
L (^> £o)> we assume a periodic lattice of such wells. A 
two-dimensional sketch is shown in Fig. |l|c. We define 
the occurrence of a nucleation event as the escape of the 
particle from the region r < s. Initial recombination 
denotes recombination with the original well, while we 
define here homogeneous recombination as capture by a 
neighboring well (Fig. |f|d). We emphasize that homoge- 
neous recombination and initial recombination differen- 
tiate histories of particles with separations larger than s. 
A homogeneous recombination event is thus represented 
by a trajectory that starts from s and ends up at a differ- 
ent well, whereas initial recombination is represented by 
a trajectory that starts from s and after excursion in the 
flat potential region returns to the initial well. Initial re- 
combination refers thus like homogenous recombination 
to particles that have completely dissociated. The only 
difference is their different subsequent history due to dif- 
fusive motion. 

In order to determine the nucleation rate, we consider 
first a single well. The bare potential U(r) does not ex- 
hibit a saddle and the particle feels an attracting force 
which, however, becomes (exponentially) small with in- 
creasing distance from the well. Once the particle is far 
enough away, it behaves thus purely diffusive and must be 
considered to be free. As in standard nucleation theory, 
the solution of the Kramers problem must obey absorbing 
boundary conditions outside the well. The introduction 
of the point r = s on the flat region where these bound- 
ary conditions are applied defines the size of a fictitious 
saddle separating the bound state from the free states. 
The exact location of this point is not determined due 
to local translational invariance in the region of constant 
potential. However, besides satisfying L ^> s S> £0, it 
should be fixed by physical conditions. A typical exam- 
ple of an experimental condition which provides a certain 
value of s is the finite resolution of an instrument which 
counts the free particles ||] . 
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To calculate the nucleation rate, we solve the station- 
ary Smoluchowski Eq. (|l|) in radial coordinates. Due 
to rotational symmetry there is no dependence on an- 
gle coordinates for d > 1. With an ansatz P(r) — 
j3{r) exp(—U(r)/kT) for the probability density, one finds 
from Eq. ([!]) for the radial current density in d dimen- 
sions 



j r = -Dexp(-U(r)/kT)d r /3(r) 



C 

„d-l 



(6) 



Here, C has to be determined and is related to the nu- 
cleation rate by J nuc = OrfC with a\ = 2, 0,2 = 2ir, and 
03 = 47r. The density is assumed to be normalized 



P(r) a d r d - L dr = 1 



(7) 



and to satisfy the linear, homogeneous, and mixed 



boundary conditions at s |22 



API 



PI, = 



(8) 



Some comments concerning this boundary condition are 
in order (see, e.g., ||). First, x m = cty/D/j is the 
Milncs length which is a measure for the nonequilib- 
rium boundary layer at an absorbing boundary. Indeed, 
•Em OC ^kT/m/j can be interpreted as the mean free 
path of a particle with thermal velocity. In one dimen- 
sion, a — £(—1/2) w 1.46. In higher dimensions, a is 
changed (e.g., due to the finite curvature of the bound- 
ary), but here we are not interested in details concerning 
this problem and we will use a as a given parameter. In 
the limit 7 — > 00, an absorbing boundary at s implies 
P(s) = 0. On a microscopic length scale, the bound- 
ary condition (||) describes the solution asymptotically 
far away from the absorbing boundary, r < s — x rn . The 
boundary layer itself (s — x m < r < s) shows a more 
complicated structure p2^p3f . Hence, Eq. (||) cannot be 
used if x m is of the order of s, or if the boundary layer 
leaks into the well. The mixed boundary condition (||) 
serves only to investigate how the results presented below 
change as we depart from the regime of strongly damped 
motion; it should be clear that for weak friction, it is not 
sufficient to consider the Smoluchowski equation but in- 
stead one has rather to use the full Kramers equation. 

Let us now continue with the derivation of the nucle- 
ation rate. Integration of (^J) together with ((7J) and (||) 
yields 



d > 1, the rate cannot be determined by using the volume 
of the symmetry group in the order-parameter (r) space 
p4j . This works only for a saddle with a well defined 
sharp maximum. We anticipate that the nucleation rate 
(0) depends in general on s. The discussion of the de- 
pendence of the nucleation rate on this parameter is our 
main goal. In the sequel the high and low temperature 
limits will be discussed separately. 



A. High-temperature limit 

For very high temperatures kT 3> Eo , the particle sees 
an overall flat potential and does not feel the well. One 
finds from M) a rate 



Jnuc o 



2d 



1 + 2{x m /s) 



(10) 



„d-l 



dre -U{r)/kT 



-,1-d 



-U(s)/kT 



,1-d 



dz 



In the high-friction limit this describes just the spatial 
diffusion out of a region r < s. In the low- friction limit, 
on the other hand, the rate is dD/(sx rn ) — \/kT /m/ see. 
Although this result cannot be right as discussed above, 
let us compare it with an estimate of the zero-friction 
result obtained directly from phase-space considerations. 
In the limit x m ^> s, the probability density in phase 
space, f(r,v), must be considered rather than the con- 
figuration space density P(r). For a constant poten- 
tial, we take / = Aexp(—mv 2 /2kT) for radial velocities 
pointing away from the well, and / = for velocities 
pointing towards the well. In the following, we restrict 
ourselves to the one-dimensional case. The density in 
configuration space (r > 0) and the current density are 
P = Jdvf = Ay/nkT/2m and j = J vdvf = kTA/m, 
respectively. Normalization implies P = l/2s and we 
find J nuc = 2j = (^/2/^)( v / kT/m)/s. Up to a constant 
factor of order one this is in rather good agreement with 
the result derived from Eq.(pTj|). 

Let us briefly determine the recombination time t± in 
the overdamped one-dimensional case. Due to equipar- 
tition at high temperatures, the ratio of free particles 
(r > s) to particles which are not free (r < s) is given 
by Nf Tee /N wc \\ — (L — s)/s « L/s. Now, the total nucle- 
ation rate and the recombination rate balance each other, 
iVwciiiuc = Nfr e<i /n. Using Eq. (JToJ) , this immediately 
yields the recombination time t\ = Ls/2D. In part C 
of this section, we will discuss the recombination times 
more deeply. 
1 



e -U(z)/kT 



(9) 



Alternatively, we could derive this result by evaluating 
the mean first passage time (see, e.g., jll|]) of a particle 
starting at r = and reaching r — s, with appropriate 
boundary conditions. 

Despite rotational symmetry of the saddle manifold for 



We notice that in the high temperature limit the as- 
sumption of a diluted gas of nucleated structures, i.e. 
the assumption L > ^0, breaks down in many physical 
applications. This is because L~ d is often a density of 
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nucleated structures which is itself proportional to an Ar- 
rhenius factor, and becomes large at high temperatures. 
In that case, initial recombination and homogenous re- 
combination can no longer be clearly separated. 



B. Low-temperature limit 

In this case, the d-dimensional probability density in 
the well is a normalized (see Eq. (|7j)) equilibrium distri- 
bution 



P (r)=A d exp{-U(r)/kT) 
with the normalization constant 



A,= 



2 \ d / 2 



exp(-E /kT) 



2nkT J 

The low-temperature nucleation rate is given by 



Jnuc — a d,AdD 



r^dr + x m s 1 - d 



(11) 



(12) 



(13) 



where £o is the size of the well. It is interesting, that 
in one space dimension (d = 1), the result ( |l3|) can be 
written in the form 



/2vr 



a s - £ 



-Jtsi 



(14) 



where Jtst is the transition state rate (Q). Note that 
despite of the vanishing of x m in the limit m — » 0, where 
the particle is overdamped, the rate remains finite since 
Jtst diverges. In the overdamped limit and for £o -C s, 
the rate is proporional to 1/s and can be expressed in 
the form: 



J =2^ 

' J nuc ^ 

7s 



kT 
2irm 



eM-Eo/kT) 



(15) 



Equation (fL5|) is the central result of this work. Accord- 
ing to Eq. (|15|), the characteristic decay time scales lin- 
early with s. In Sect. IV we show that the equilibrium 



nucleation reate of kinks has the same dependence on s. 

For weak damping x m 3> s, we find that J nuc devi- 
ates from the transition state result only by a factor of 
1.73. This (relatively small) deviation is due to the above 
mentioned fact that the inner structure of the boundary 
layer must be taken into account and cannot be described 
simply by a boundary condition. One finds exactly the 
transition state result by an appropriate matching of the 
probability density P = AyJirkT /2m (see end of the last 
subsection) to the part of the solution ([ll]) associated 
with right-moving particles (r > 0). Hence, our result 
has a behavior similar to Kramers' moderate-to-strong 
damping formula for the quadratic saddle. As mentioned, 
( |ll| ) is a wrong distribution at low damping since local 



equilibrium in the well is not established, and the nucle- 
ation rate at 7 — > is given by Eq. (||). 

In two and three space dimensions, the rates can be 
expressed in the form 



7 = 



and 



' J nuc ~ 



^0 



exp(-£ /fcr) 
7 ln(s/£o) + Xm/s 



2fo 



7 V 2nkT 1 - (&/*) + (CoW* 2 ) 



(16) 



eM-Eo/kT) , 



(17) 



respectively. In contrast to the one-dimensional case, the 
well size £0 cannot be neglected for d > 1. Moreover, in 
two dimensions the s-dependence is logarithmically weak, 
and in three dimensions its influence can be neglected. 
We mention that a discussion of the result in the limit 
£o^0 requires some care, since £o> w o> an d Eq are usu- 
ally not fully independent. 

Furthermore, in d = 2, 3, neither £0 nor s are 
necessarily equal to the usual transition state ob- 
tained by extremalization of the effective potential 
kT\n{r 1 - d e^{U/kT)) = U(r) - kT(d- 1) ln(r). In con- 
trast to the usual quadratic saddle, it is not only the 
vicinity of this transition state which contributes to the 
integral in Eq. (^|) . It is a larger region around the tran- 
sition state that counts, of which £0 an d s, respectively, 
are the lower and the upper boundaries of inegration. 



C. Lifetime of free particles 

We return now to the periodic lattice of the wells 
sketched in Fig.|l|c. Each well can then be associated 
with a unit cell of volume Vg — L d much larger than 
the volume associated with bound particles (oc s d ). At 
equilibrium and at low temperatures, kT <C Eq, the ra- 
tio of particles in the well to free particles is given by 
N we i\/Nf Tee w 1/VdAd- In steady state, the generation 
rate N wc i\J nuc equals the homogeneous recombination 
rate iVf ree /r. The lifetime of free particles is thus given 
by 



Td 



V d A d 

J n.ii.a 



(18) 



For very strong damping and for £0 "C s, Eqs. (|1J)-([T^ 
imply 



n 



Ls 
2D 



To = In — 

2nD VCo 



T3 = 



(19) 
(20) 
(21) 
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In one dimension, the lifetime scales linearly with s, while 
the s dependence in higher dimensions is weak or negli- 
gible. The higher the dimension the more unlikely is 
the probability to reach the original well by pure diffu- 
sion. This is related to a divergence of the lifetimes as 
£o vanishes. We also mention that ( pll) corresponds to 
the standard result in the theory of diffusion controlled 
reactions, which states that the rate per volume is given 
by 47r£>& (see, e.g., §). 

For the one-dimensional case, the lifetime can be cal- 
culated directly by solving the Smoluchowski equation 
(|j) for a flat potential, for absorbing boundary condi- 
tions -P(O) = P(L) = 0, and for an injected current 
j(s + 0) — j(s — 0) = j nuc at r = s. If P(r) is normal- 
ized, the lifetime is then given by 1/ 'jnuc- One finds that 
recombination can be understood as a sum of two contri- 
butions, i.e. it holds j nuc = j s + jh, where jj 1 = Ls/2D 
associated with (|l|) and j^ 1 = L(L — s)/2D associated 
with diffusion to the neighbor well. 

Below, in the discussion of equilibrium kink nucleation, 
we will first calculate the lifetime independently from the 
nucleation rate and from equilibrium statistical mechan- 
ics. Using then both results for the nucleation rate and 
the lifetime, we will derive the equilibrium kink density 
and show agreement with the statistical mechanics result 
for the kink density. 



D. Effect of a weak force 

In this section we discuss the influence of an external 
weak force F on nucleation in one dimension and at low 
temperature. The potential U (r) is now replaced by the 
new potential U(r) — Fr. We are interested in the re- 
sponse to leading order with respect to F. The saddle is 
still flat enough such that diffusion dominates the drift in 
the relevant region: kT 3> Fs. For finite F, the current 
to the left (j-) and to the right are no longer equal 
and the nucleation rate is J nuc = j+ + j- . Interestingly, 
as we will show below, the transport current Jtrans is de- 
termined by the imbalance of the forward nucleation j + 
and the backward nucleation j_ . To leading order in F, 
the transport current is thus given by 



3+ - 3- 



(22) 



The symmetry relation U(r) — U(—r) implies j+(F) — 
j-(—F), and J nuc is an even function of F while Jtrans 
is an odd function of F. We find 



.]+ 



D 



A X F 
~~kT~ 



1 - exp(-Fs/kT))~ 



(23) 



where Ai is the normalization constant given by Eq. fll2| ) 
with a renormalized frequency and an energy of the new 
minimum 



U 




2to 3 w§ ' 



(24) 



Eo = Eq I 1 



F' 2 



(25) 



In Eq. (||) we defined = d$U\ Q . The nucleation 
current J nU c and the total current Jtrans can then be 
written in the form 



J n 



FAi 

TO7 



c ° th G§0 

FAi 



J i, 



(26) 
(27) 



respectively. This yields a relation between the two cur- 
rents, Jtrans = Jnuc ta,vh(F s / 2kT) . A similar result has 
been derived in Ref. Q for photo-generated currents. Ac- 
cording to our assumptions, the results are valid in lead- 
ing order with respect to F. An expansion gives 



J i, 



Fs 
2kf 



■J, , 



(28) 



Since the nucleation rate is proportional to s -1 , the total 
current is independent of s as it must be. 

We now present two derivations to show that the trans- 
port current Eq. ( p2| ) is indeed just determined by the 
imbalance of the forward and backward nucleation rates. 
For weak driving forces the transport current is equal to 
the density of free carriers, nf ree , multiplied by their drift 
velocity u = fiF, 



U71.fr, 



/LtFrifre 



(29) 



As usual, the mobility and the diffusion constant are 
related via the Einstein relation D = fikT. The num- 
ber of free carriers is determined by the balance of the 
generation of free particles and their recapture into the 



well, J n 



A'free/r with iVf r , 



Here, J n 



determines the frequency with which free particles are 
generated, and iV wo n ~ 1. The lifetime t which de- 
termines the recapture back into the well is a function 
of F. But to leading order in F, the equilibrium life- 
time Eq. (|I<| ) is all we need. Eliminating the density 
of free carriers and using the Einstein relation, we find, 
Jtrans = ^FrJnuJL = (Fs/2kT) J nuc . This is in accor- 
dance with Eq. (p8|), and thus Jtrans = j+ — j- holds. 

A second derivation of this result proceeds as follows. 
We still consider a periodic potential with period L and 
which is symmetric around the origin. It is sufficient to 
consider the range < r < L. We view the transport 
current as a consequence of the source currents j+ at 
r = s and j- at r = L — s. The wells act as particle ab- 
sorbers. The solution to this problem is a superposition 
of solutions to two problems each with one source alone. 
The current which flows back from r = s in the pres- 
ence of the source j + is denoted by j+ and the current 
which flows forward to the next well by . Similarly, 
the source j_ alone leads to a current back into the well 
at r — L denoted by j_ and a current into the well at 



G 



r = denoted by j_ . The transport current is then given 
by 

Jtrans = 3+ + 3-- (30) 
Continuity of current requires 

3 + =3+-3% 3-=3 R -J- ■ (31) 

As we have seen already the source currents have the 
symmetry j + (—F) = j_(F). The homogenous recom- 
bination currents j+ and j_ are related by symmetry 
according to j+(—F) = —j^(F). On the other hand, the 
initial recombination currents are even functions of the 
field j+(-F) = j^(F) and j+(—F) = jf(F) and more- 
over, they are equal in magnitude but differ in their sign, 

ft(F) = -j*(F) . (32) 

The initial recombination currents are maximal for F — 
0. With increasing F the recombination current jf de- 
creases because homogeneous recombination increases. 
Similarly, the initial recombination current j_(F) de- 
creases because there are fewer carriers activated into the 
high energy region of the potential. Using Eq. ( pCf ) and 
Eq. @ we obtain J trans = j+-j-+j^{F)+j R {F). But 
as we have seen the sum of the two initial recombination 
currents cancel one another and thus Jtrans = j+ — 3- ■ 
Thus at low fields the difference of the two activation 
rates directly determines the transport current. This can 
also be shown by a direct calculation. 

The direct relationship of the rates j+ and j_ to the 
transport current demonstrates that these are physically 
meaningful and useful quantities. 

IV. EQUILIBRIUM KINK NUCLEATION 

In this section we investigate the dynamics of a string 
of particles coupled to each other harmonically and mov- 
ing in a sinusoidal potential. The particles are subject 
to damping and noise and might be subject to an ex- 
ternal driving force. This model is known as the driven 
and damped sine-Gordon chain. It has a long history 
and due to its wide range of applicability, from kinks in 
surface steps on various materials to the motion of flux- 
ons in long Josephson junctions, has been widely stud- 
ied [ pl]Jl^j25| -p5[ . In the overdamped limit, which is of 
interest here, there are only two types of elementary ex- 
citations. There are small amplitude phonon like exci- 
tations and, more intcrstingly, highly non-linear struc- 
tures, called kinks or solitary structures, which describe 
the transition from one valley to another. In chemical 
physics, kinks are discussed in various polymers fl36| . 
Our concern is the statistical mechanics of such a sys- 
tem which we take to be so large that one can define a 
density of kinks and anti-kinks. Of interest is a theory of 
the thermal equilibrium nucleation of kink-antikink pairs 
and particularly the role of initial recombination f|] • We 



show that a theory of equilibrium kink nucleation, i.e., 
for a vanishing external driving force, must take into ac- 
count initial recombination. 

At small temperatures the kink-antikink gas is diluted, 
and the kink density sets an upper length scale over which 
nucleation and annihilation processes have to occur. In 
the framework of equilibrium statistical mechanics, kinks 
and antikinks are regarded as free particles. Below we 
shall explain that at equilibrium the distance between 
a 'mathematically unbounded' kink-antikink pair is in- 
finitely large. Apparently, this seems to contradict the 
notion of a free kink in a finite system. For reasons of 
consistency, it is thus necessary to develop a picture of 
the nucleation and annihilation processes which permits 
essentially free diffusive motion during the lifetime of a 
kink. 

The nucleation, the dynamics, and the recombination 
of kinks and antikinks in space-time is schematically il- 
lustrated in Fig. 1 of Ref. jjj . Out of equilibrium, in the 
presence of a strong force (see Fig. la of Ref. |ffl]), kink 
and antikink are driven apart after a nucleation process 
(empty triangles) by the force and annihilate eventually 
with an antikink and a kink originating from a differ- 
ent nucleation process (rectangles). This picture has to 
be qualitatively modified in the equilibrium case without 
force. As shown in Fig. lb of Ref. jy , the diffusive motion 
of the free kinks gives them a strongly enhanced probabil- 
ity of initial recombination (closed loops). Neglecting the 
closed trajectories is inconsistent with the experimental 
definition of free kinks. Former works f^,0 (see also Ref. 
]j"l|) do not describe the initial recombination processes 
and count only the negligibly small fraction of extended 
trajectories (see Fig. lb of Ref. jlj]) which describe homo- 
geneous recombination. Consequently, these works arrive 
at a much too low nucleation rate with an activation en- 
ergy associated with a kink triple, instead of the 
pair energy 2Ek- Hence, these works predict an incom- 
prehensible breaking of the kink-antikink symmetry. A 
further consequence is a mean kink lifetime proportional 
to exp(2£ fc /fcT) @||,||, whereas our theory leads to a 
kink lifetime proportional to exp(EkfkT). 

To be specific, let us consider now the overdamped 
sine-Gordon equation j?7| 

7 «9 t = -V o sin0 + F + Kd%6 + ( , (33) 

which describes the dynamics of an order-parameter field 
(a string of particles) 6(x, t) in a periodic potential of am- 
plitude Vq, and with a coupling constant n. Unless other- 
wise stated, the force F is set to zero. We assume periodic 
boundary conditions 9{C + x,t) = 8(x, t), where C is the 
sample length which exceeds every other relevant length 
scale of the problem (except the diverging size of the 
mathematical critical nucleus). Furthermore, C, denotes 
a weak white noise force associated with the temperature 
T, i.e. with zero mean (0=0 and the correlation func- 
tion (C(x,t)C(x,t)) = 2jkT5(x-x)8(t-t). The uniform, 
stationary, and linearly stable states are given by 9 St i = 
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2lw (Peierls valleys) with integer I, and have equal ener- 
gies. There exists an energy functional E[6] such that Eq. 
( |33| ) can be rewritten in the form jdtO = —6E[6]/66. For 
a weak finite force, two adjacent Peierls valleys are sep- 
arated in function space by a saddle which corresponds 
to a kink-antikink pair. A kink 9 k (x — xq) centered at Xq 
connects a Peierls valley 9 s ,i with its neighbor 9 s ,i+x- An 
antikink is reversely defined by 9 a = 9 k (—x+xo). A kink- 
antikink pair at location Xq and with a (not too small) 
separation r can be written approximately as 9n(x) = 
6 k (x-x +r/2)+9 k (~x+x +r/2)-2TT(l+l). In the pres- 
ence of a weak force F, the mathematically exact saddle 
point of the energy functional corresponds to a pair with 
a separation £ = £ ln(Vo/.F) where £o = y/n/Vo is now 
the kink size. In the equilibrium limit where F — ► 0, the 
separation £ exceeds the inverse kink density n _1 or even 
the system length C. In this case, the mathematically ex- 
act critical nucleus has no physical meaning. In order to 
overcome the problem of the large critical nucleus, we 
introduce an effective critical nucleus with a separation s 
of the kink and the antikink which is larger than a kink 
width but much smaller than an average equilibrium sep- 



aration L 



of kinks. The effective critical nucleus 



corresponds to the flat barrier with size s, in analogy to 
the effective saddle discussed in the previous sections. 



variables 770 and 771 are the orthonormal-mode coordi- 
nates which belong to the kink-separation mode and to 
the translational mode, respectively. It holds 



drfc = dr 2 / (59 N /Sr) 2 dx 



di] 2 = dxl 



In /5xo) 2 dx 
(35) 



Using the quasi-Goldstone and the Goldstone modes 
given above, one finds for the ratio of the integrals in 
Eq. (|34| ) a value 2L/s. The normalized partition func- 
tion of the damped degrees of freedom at the saddle, 



Zn 
Zm 



I 

2^ 



\ 



A A x 



n 



(36) 



contains the stability eigenvalues X^' N of the metastable 
state (index M) and the critical nucleus (index N) with 
respect to perturbations with a decay (or a growth) oc 
exp(Ai). The (quasi-) zero modes are excluded in the 
products. For a well separated pair, Eq. ( |36| ) is the 
renormalized partition function of a kink-antikink pair 
without self-interaction and is given by the square of the 
single kink partition function. Hence, Zn/Zm = 4T/2n, 
where r = V0/7 J30| . The kink nucleation rate per length 
finally becomes 



A. The nucleation rate 

We derive now the kink nucleation rate per length, 
Jnuc- Note, that due to the periodic array of Peierls val- 
leys, a kink-antikink and an antikink-kink are equivalent 
excitations. If the nucleation rate of the former and the 
latter are denoted by j+ and j_, respectively, the rate 

becomes J nuc = j+ +j Of course, at equilibrium it 

holds j+ = j-. The presence of spatial degrees of free- 
dom requires a multi-dimensional Kramers theory p4[ . 
The relevant coordinate is the kink-antikink separation 
r which can be associated with a quasi-Goldstone mode, 
besides the usual translational Goldstone mode of the 
pair. The quasi-Goldstone mode is 59jy/5r = 9' k (x — xo + 
r/2)/2 + 0' k (—x + xo + r/2)/2 associated with an infinites- 
imal variation of r. The (orthogonal) Goldstone mode 
associated with an infinitesimal displacement Sxq of the 
pair is S9n/Sx = 9' k (—x + x + a/2) - 9' k (x — x + s/2). 

The stationary Fokker-Planck equation associated 
with the Langevin equation (|33| ) is solved in an analo- 
gous way as m Sect. |n| §. Dc grees of freedom which 
are transverse to the two Goldstone modes can be inte- 
grated out by Gaussian integration |3C[ ], The rate per 
length of a kink-antikink pair becomes then [EJ 



= }_Zn Jo 
J+ L Z 



m(L) 



dru 



Jo dr lo 



exp(-2E k /kT) 



(34) 



One concludes an activation energy Eo — 2E kl where 
E k = 8\/kVq is the equilibrium- kink energy [BO . The 



4r 

Jnuc = 2j+ = 2— exp(~2E k /kT) 

ITS 



(37) 



As discussed in Sect. (Ill), the 1/s-dependence indicates 



the strong dependence of the rate on initial recombina- 
tion events in one-dimensional systems. 



B. Lifetime of kinks 

In the following, we show that the result ([37]) is con- 
sistent with equilibrium statistical mechanics, which pre- 
dicts an equilibrium kink density p0[ 



2V E k 



eM-Ek/kT) 



(38) 



In order to derive this result independently, we will cal- 
culate first the kink lifetime r. By using the balance 
equation J„ uc t = n we obtain then a kink density which 
has to be compared with (38). 



The kink lifetime r for fixed antikink density n can be 
calculated with the help of a Langevin equation for the 
kink separation r. This Langevin equation follows from 
a projection of the sine-Gordon Eq. ( |33"1) onto the quasi- 
Goldstone mode discussed above. The Fokker-Planck 
equation which is equivalent to the Langevin equation, 
reads in the stationary case d r (FP — Dd r P) = 0. Here, 
we assumed a finite force F for later use. The effec- 
tive force acting on the separation coordinate is given by 
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F = 2(iF with a mobility fi = 2'KK/"fEj Cl and the diffusion 
constant is given by D — 2/ikT/2iT |2^,|3(|. The values of 
F and D for the relative coordinate r are twice as large as 
for a single kink. Now, we return to the case F = 0. The 
stationary Fokker-Planck equation must be solved with 
a source at r — s and with sinks at r — and r = n^ 1 . 
The source describes the nucleation of a pair, and the 
sinks model kink-antikink annihilation. The mean kink 
distance n~ l corresponds here to the well distance L in- 
troduced in Fig. 0c for one dimension. 



Sect. Ill D the transport current can be directly related 



Now we can proceed as in subsection [II C. Integra 



tion of the stationary Fokker-Planck equation leads to 
a piecewise constant current density, whereas the source 
implies a discontinuity of the current density of strength 
j + at s. The absorbing boundary conditions demand 
P(0) = P(n- r ) = 0. The lifetime r is defined by the 
ratio of the total probability J dsP and the injected cur- 
rent j + . One finds then immediately the result ( |l9| ) with 
L replaced by 1/n, i.e. 

(39) 



2Dn 



The balance equation together with ( p7| ) implies imme- 
diately the equilibrium density (|3^) which proves consis- 
tency with statistical mechanics, i.e., n = n eq . 

As it must be, the stationary kink density ( |3S| ) is inde- 
pendent of the specific value of s. The s dependence of 
the rate and of the lifetime can be illustrated in Fig. 1 of 
Rcf. 0| . The kinks which are to be counted in [0, C] at a 
fixed time t must have been nucleated in a strip (t — t, t) 
of width r oc s. Since the number of counted kinks has 
to be independent of s, the density of generation events 
in this strip must be proportional to 1/ s. A variation of 
s affects only kinetic quantities like the time scale, but 
not thermodynamic equilibrium quantities like the kink 
density. 



C. Effect of a weak force 

In the presence of a weak force F, we expect a drift in 
the order parameter field 9. Kinks and antikinks which 
pass a given location with average velocity u and —u, 
respectively, each advance the field by 2tt. Thus the av- 
erage speed of the order parameter field is 



(d t 0) = 2nu(n + m) 



(40) 



where m and n are the average kink and antikink densi- 
ties. In our problem it holds n = m. To leading order in 
the field F the velocity is simply determined by the kink 
mobility u = /j,F and the kink density is determined by 
the equilibrium density. Thus to linear order in the force 
F the order parameter has a velocity flol 



(d t 6) = 47TfiFn eq 



(41) 



We will re-derive this result directly from the kink nu- 
cleation rates. Like in the single-particle problem fl2q) of 



to the asymmetry of the kink nucleation rates. The kink 
nucleation rate J nuc is balanced by kink recombination, 
Jnuc — n/r. Here t is the lifetime of kinks, which as we 
have seen is at low fields determined predominately by 
initial recombination and at F = is given by Eq. (|3S|). 
Eliminating the kink density from Eq. (f40|), we obtain 
(dtO) = AiruTJ nuc . Now we find that the difference of the 
two nucleation currents is given by 



3+ - 3- 



2D n 



(42) 



To leading order in the field, it follows immediately that 
the average velocity of the displacement field is directly 
determined by the imbalance of the nucleation rates, 



(d t e) = —( J+ -3-) ■ 

This expression is equivalent to Eq. (Bl 



(43) 



V. CONCLUSION 

In this paper we have shown that diffusion controlled 
initial recombination strongly affects the activation rate 
in problems with flat potentials. The effect is partic- 
ularly important in one-dimensional systems where the 
probability of diffusion back to the original well is very 
large. The absence of a well-defined barrier maximum 
requires the definition of an effective size of the barrier, 
i.e., a location where for the evaluation of the rate ab- 
sorbing boundary conditions to the probability density 
are applied. In one-dimensional systems the time and, as 
a consequence, the inverse rate, scale linearly with this 
size. 

Furthermore, we discussed equilibrium kink nucleation 
as an example where diffusion controlled initial recombi- 
nation determines the kinetics. The effective saddle cor- 
responds here to an effective critical nucleus. A certain 
arbitraryness of the nucleus size reflects the uncertainty 
of the notion of a free kink. A specific choice of this 
size, however, has to be related to physical considera- 
tions: obviously, it has to be larger than the kink width 
but much smaller than the mean kink distance. We have 
shown consistency of our results with equilibrium sta- 
tistical mechanics. Furthermore, we have shown that the 
activation rates directly determine the transport current. 
Our work not only demonstrates that rates and lifetimes 
are useful quantities even if they depend explicitly on 
a length scale which separates free particles from bound 
particles. It also demonstrates that an evaluation of such 
rates is necessary to provide a physically meaningful dis- 
cussion of problems in which diffusion controlled initial 
recombination plays a dominant role. 
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FIGURE CAPTIONS 



FIG. 1. a) Potential well with an activation barrier at r = s 
and with an activation energy Eo; b) Potential well without 
a barrier maximum, but with a flat region; this case has to 
be treated by introducing a threshold point s. c) Lattice of 
potential wells with distance L. d) Typical trace of initial 
recombination (ir) and of homogeneous recombination (hr) 
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